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Abstract 

We investigate using plane fronted gravitational wave space-times as model systems to study 
loop quantization techniques and dispersion relations. In this classical analysis we start with planar 
symmetric space-times in the real connection formulation. We reduce via Dirac constraint analysis 
to a final form with one canonical pair and one constraint, equivalent to the metric and Einstein 
equations of plane-fronted-with-parallel-rays (pp) waves. Due to the symmetries and use of special 
coordinates general covariance is broken. However, this allows us to simply express the constraints 
of the consistent system. A recursive construction of Dirac brackets results in non-local brackets, 
analogous to those of self-dual fields, for the triad variables. Not surprisingly, this classical analysis 
produces no evidence for dispersion, i. e. a variable propagation speed of gravitational pp-waves. 



PACS 04.20.Fy, 04.30.Nk, 04.60.Ds, 04.60.Pp 



I. INTRODUCTION 

The success in constraining modified dispersion relations [l-4] has renewed efforts to see 
whether, in the context of various approaches to quantum gravity, such modifications arise. 
This is interesting even in model systems where the quantization may be more unambigu- 
ously carried out and where it is possible to identify the origin of the modifications, should 
they appear. For instance, this has been explored in the context of polymer quantization of 
scalar fields in flat space-time {5, 6]. In this case the origin of the modifications lies in the 
choice of classical polymer variables, in particular the length scale required to express the 
exponentiated momentum variable, rather than in a granularity of spatial geometry. While 
loop quantum gravity there are heuristic results suggesting that there might be modifications 
to dispersion relations |7-9|, it would be interesting to investigate possible modifications in 
a model system in which both the origin of the modification is clear and in which the quan- 
tization may be completed. This paper explores whether the symmetry reduced space-times 
of plane-fronted gravitational waves with parallel rays may be a suitable context in which 
to explore modifications to dispersion relations. 

Classical plane gravitational waves are, like homogeneous and isotropic cosmological mod- 
els, among the most simple exact solutions of General Relativity (GR). Despite the nonline- 
arity of GR, due to the symmetry of the model pulses of plane fronted with parallel rays (pp) 
waves travel without dispersion and leave space flat outside the pulse; they form a "wave 
sandwich" with the gravitational wave pulse between regions of flat space. In fact the only 



gravitational waves with this flat space sandwich property are pp- waves [10] . Due to their 
simplicity pp-waves promise to be good candidates to test quantization techniques for pure 
gravity. It is also intriguing that the Einstein equations for plane electromagnetic waves 
coupled to gravity take the same form as the pp-waves 11], suggesting that a quantization 



and study of dispersion relations of pp-waves could be extended to this Einstein-Maxwell 
theory model. So the quantization of pp-waves may yield an answer to whether quantum ef- 
fects of 3-dimensional geometry can lead to dispersion of gravitational (or electromagnetic) 
plane waves, whether spatial granularity can yield an energy-dependence of the speed of 
gravitation waves or light. 
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In view of this eventual goal, we formulate polarized parallel plane waves in terms of the 
real connection variables, proceeding from rather general assumptions about homogeneity 
in two dimensions, to a form equivalent to a standard form in the literature, given below. 
Despite the simplicity of this well-known result, the canonical way to this goal is not trivial. 

The metric of pp-waves propagating in z direction, in a "Rosen-type" chart, given by 



Misner, Thorne, and Wheeler 



111 is 



ds 2 = _ dt 2 + L 2 e W dx 2 + L 2 e -W dy 2 + dz \ (1) 



(See also Ehlers and Kundt [12|].) This metric has a convenient interpretation for our pur- 
poses. The function L, called the "background factor" , is determined by the free function j3, 
called the "wave factor". Both L and /3 are functions of v := t + z or u := t — z. L satisfies 
the Einstein equation 

L" + (f3') 2 L = 0. (2) 

(In this equation the prime denotes a derivative with respect to u or v.) This single equation 
survives the reduction of GR. In it 0' acts as a "time-dependent" angular frequency. In the 
light cone coordinates u - once (3(u), L(0) and L'(0) are specified - the function L{u) is 
determined just as in a simple 1-dimensional mechanical system. 

The above form of the metric as well as equation ([2]) are valid only for waves traveling in 
the positive or in the negative z direction. The combination of both, i. e. colliding waves, 
are a problem of a much higher degree of difficulty jl^, [l^] . 

An important feature of the above coordinate system is that it does not globally cover the 
space-time inhabited by plane waves. One may see this even with a short pulse, with /3 ^ 
only in an interval (u_,u + ). For example choosing L = 1 in the region of flat space in front 
of the pulse, at the location of the pulse, where (3 ^ 0, L" becomes negative due to fl2]) and 
so L decreases inside the wave, as long as L > 0. (We assume the pulse to be short enough 
and not too strong so that L > everywhere inside. For details of this approximation see 



11].) In the flat space behind the pulse, when (3 is constant again, L" = and L is a 



linear function, which has to be matched smoothly to L at the end of the pulse. This leads 
necessarily to L = at a certain value u 2 > u + somewhere behind the pulse, a coordinate 



singularity in flat space jll[ 
u < u 2 - 



3. 



In this case the coordinate system is valid in the region 
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One can choose - and we shall do so - L to be a non-constant linear function on both 
sides of the wave pulse, with one zero in front, at u — Ui, U\ < U-, and one behind it. In 
this case, the coordinates cover a slice of the gravitational wave sandwich. In detail, we have 
the three sub-intervals of the coordinate range (wi,^) 

(1) Ui < u < U-, V — const. > 0, f3 = 0, flat space in front of the pulse, 

(2) U- < u < u + , L" < 0, f3 7^ 0, the pulse, and 

(3) u + < u < u 2 , L' = const. < 0, f3 constant, flat space behind the pulse. 

The boundary conditions at u — Ui$ are flat-space boundary conditions with constant (3 
and L = at the coordinate singularities. Despite the regions of flat space before and after 
the pulse, neighboring test particles in the xy plane accelerate and fall towards each other 
as the wave passes |ll|, 



1Q. 

In the next section we review the canonical variables and the polarized Gowdy model of 
[if) ] . The consistent reduction to the pp-wave case is accomplished in sections III Dl and III El 
The Dirac brackets are constructed in section HTT1 Time evolution in the preferred coordinates 
is discussed in section IIII Dl The results of the classical calculations are summarized in 
section EJ Finally, a note on the orthogonality of the connection and the Immirzi parameter 
is inHVl 

II. SYMMETRY REDUCTION 



A. The connection variables 

We formulate the system, after a 3+1 split, with the usual densitized triads E a i and con- 
nection components A a % of the real connection formulation. (See Ref. {^(j] for a review.) We 
denote spatial and % = 1, 2, 3 as an su(2) Lie algebra index. By homogeneity, 

we further assume that on every spatial slice they are functions of z alone, with their time 
dependence to be determined by equations of motion. The system has a close formal analogy 



to the polarized Gowdy T 3 model analyzed by Banerjee and Date 16|, [17|- For this reason 



we have chosen essentially the same notation, so that at many points we can refer to these 



papers. Of course, the angular variable 9 of 16( had to be changed to z and, as is clear from 



this notation, we are also not working with a compact spatial topology. 
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The symmetry reduction from full GR to the Gowdy model is carried out in a process 



outlined, for example, in the appendix of 



18]. In the present case the symmetry reduction 



is briefly the following: We start with a principal fibre bundle of SU(2), the gauge group 
of loop quantum gravity (LQG), over a three-dimensional space manifold. The spatial 
symmetry in the presence of pp-waves consists of translations in the x and y direction, 
the orbits of which are planes parallel to the xy-plane. The space manifold is decomposed 
into an orbit bundle with a one-dimensional basis manifold, the z-axis, called the reduced 
manifold. The "reduced bundle" is the trivial principal SU(2) bundle over a single coordinate 
neighborhood z\ < z < z 2 of the reduced manifold, where Z\ and z 2 correspond to the null 
coordinate boundaries Ui and u 2 in the introduction. The symmetry reduction proceeds 
with a decomposition of the bundle connection. The latter one is separated into the reduced 
connection, i. e. the restriction along z, A z \z)ri, and into scalar fields on the reduced 
manifold, A x l (z)Ti, A y \z)Ti where Tj are the usual su(2) generators. By a choice of gauge 
the reduced connection is assumed to lie in the subalgebra generated by T3, the scalar fields 
in the subspace spanned by T\ and t 2 , so that the matrices A a l and the canonically conjugate 
densitized triad matrices become block-diagonal [19] 

= E? = 0, Aj = A p 3 = 0, with p = x,y, I = 1,2. (3) 

As usual in LQG, the connection A& is defined as the combination 

a: = y:- 1 k: (4) 

of the torsion-free spin connection IV and the extrinsic curvature K a l . The Barbero-Immirzi 
parameter is denoted with 7. These variables are subject to the usual constraints of canonical 
GR, the Gaufi, the diffeomorphism, and the Hamiltonian constraint. 

Due to the planar symmetry of the waves, the phase space variables are free of dependence 
on x or y. We restrict integration in the xy plane to a finite, fiducial patch with area A a . 
Integrating the symplectic structure over this patch gives 

Q = — fdz (dA 3 z A d£g + dA? A dE p j) (5) 

where k' = AttG. For the rest of this article we use k = k'/A . Following [16| we denote E z 3 
by £ and A 3 by and introduce polar coordinates in the "1-2" plane, 

E x 1 = E x cos (3, E X 2 = E x sin(3, 

(6) 

E v i = -Ev sin p, Ey 2 = E^ cos /3, 
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A x = Aj;COs(a + /3), A 2 = A x sin(a + f3), 

(7) 

Ay 1 = —A y sin(a + (3), A y 2 = A y cos(a + (3). 

The canonically conjugate connection variables to the radial variables E x and E y , [3 and (3 
are 

1 1 

:= — A x cos(a) and K y := — A y cos(a) (8) 

and 

P p := -E x A x sin(a) and := A y sin(a), (9) 

respectively. 



B. The polarization condition 



Following Banerjee and Date [16|] we carry out a reduction to the polarized model by 
setting (3 — j3. This ensures that the Killing vectors d x and d y are orthogonal. This means 
that E x and E y become orthogonal in the sense that 

E x iE\ = (10) 

and that the spatial part of the metric ([1]) becomes diagonal. Denoting the spatial distance 
as ds 2 

Tpy pi px rpy 

ds 2 = £±L dx 2 + S ±_ dy 2 + f_±L d ^2 (n) 

E* E» E 

At this point we do not yet specify the lapse function and the shift vector. After redefinition 
of the angular variables and their momenta 

p/3 _ pP 

£:=P-h P*-= g (12) 

and 

77 := /5 -h /5, P":= I (13) 

the condition (3 = (3 may be imposed in the form of a second-class constraint, ^ := /3 — /3 = 0, 
whose Poisson bracket with the Hamiltonian constraint does not vanish weakly. When {£, H} 
is added to the constraints, then together with £ it forms a pair of second-class constraints, 
weakly Poisson-commuting with all the other constraints. After introducing Dirac brackets, 
this pair of constraints can be imposed strongly, thus eliminating £ and flo| . 



The remaining phase space variables are A, K x , K y , 77 and £, E x , E y , P v ; their Dirac 
brackets are equal to the Poisson brackets, as they all have weakly vanishing Poisson brackets 
with the two strongly imposed constraints. 



C. The gauge constraint 

We have accomplished the symmetry reduction and imposed the polarization condition. 
In these variables the GauB constraint G, the diffeomorphism constraint C and the Hamil- 
tonian constraint H reduce to (see |l6| ) 

G= — \£' + 2P% (14) 
C = - [K' X E X + K'E y -£'A+- v'P% (15) 

K y 7 



H- 1 



2k\[E 



yG 2 + (K X E X + K y E y ) (j- + 2A S j£ + 2K x E x K y E y + 



(16) 



1 I ( f /)2 _ e 2 {^r - 1 + 2 (- + |P) EPv + 8'P* + 2£P*>' 



2 r ' \E y /•; ' ) 1 { \E X E y ) 



where E := £E x E y is the determinant of the 3-metric. The prime means derivative with 
respect to z and k is the gravitational constant. The total Hamiltonian is 

H tot = A J dz (AG + nC + NH) . (17) 

The authors of [16] point out that for the polarized Gowdy model an orthogonality condition 
on the A a l analogous to fflOl) . namely 

AjAj = 7 {K x E xl - K y E yl ) = (18) 

in the above variables, is not conserved under evolution. More precisely, the Poisson bracket 
of the condition with H does not vanish weakly, it would give rise to a further constraint, 
and so on, rendering the system inconsistent. Nevertheless, a careful analysis of the spin 
connection and extrinsic curvature derived from the metric (CQ) shows that pp-waves do 
satisfy this condition. We will come back to this issue in Section IIVI 



For the present we follow the simplifications in 16] one more step by strongly imposing 
the GauB constraint together with the associated gauge fixing condition 77 = and thus 
remove the variables rj and P v . Again the Dirac brackets of the remaining variables are 
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equal to their Poisson brackets, defined by 



, . . , . '8F8G 8F8G 8F 8G SF 5G 
{ F,G} = K ldz(^--- u + — — -— — 



(19) 
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When we have carried out this, the diffeomorphism constraint fl!5p drops its last term and 
the Hamiltonian constraint becomes equal to 
1 1 



H 



2k ^E 



[2K x E x K y E y + 2(K X E X + K y E v )£A] + 



1 1 

4k y/E 



S> 



E yl E x> 

~Ev ~ ~E* 



2SS' 



, (E yl E x 



\Ev + E x 



+ (s'y + ASS" 



(20) 



This last step reduced the system to three canonical pairs, related by two first-class con- 
straints, so that the field theory has exactly one phase space degree of freedom per spatial 
point. As shown in Ref. 16] the algebra of the constraints is the correct one for canonical 
GR. 



D. Reduction to pp-waves and the spatial Einstein equation 

Now we reduce the theory to a model equivalent to the one formulated in metric variables 
by (JTJ) and 02]), again using Dirac's constraint analysis. The metric of ([T]) contains two 
functions, L and rather than the three functions S , E p in ( II ip . The one Einstein equation 
02]) is not equivalent to the remaining constraints C = or H = 0. So we need (at least) 
one more constraint. Comparing the spatial part of the metrics ([1]) with ffTT]) we see that 
we need the primary constraint g zz = 1, or 

B := S — E x E y = 0. (21) 

Of course there is no guarantee that the resulting system is consistent - after all the polarized 
Gowdy is already reduced to two phase space degrees of freedom - but this system is simple 
enough so we can introduce the appropriate constraints in the special coordinate system. 

For the reduced theory to be consistent, the new local constraint B must be preserved 
under evolution of the total Hamiltonian constraint. The Poisson bracket with the smeared- 
out diffeomorphism constraint 

C[n] := A [ dzn(z)C(z) (22) 
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IS 



{B(z), C[n}} = -n{z)B\z) + 2n\z)E x {z)E v {z). (23) 

Generally this is not weakly equal to zero, because the new constraint is not invariant under 
local diffeomorphisms due to the different nature of £ and E p . The variable £ transforms as 
a scalar, whereas E x and E y transform as scalar densities, as can be seen from the Poisson 
brackets 

{£(z), C[n}} = n(z)£'(z), {ET(z), C[n}} = (n(z)E»(z))'. 

Only when the test function n is constant along z - meaning the shift vector depends only 
on t - is the combination B of £ and E p meaningful and conserved under the action of C [n] . 

In the end the failure of B to be diffeomorphism invariant is not a surprise. Demanding 
g zz = 1 we obviously restrict local diffeomorphism invariance. (In the special case (II]) the 
shift vector is equal to zero.) The local constraint C(z), in contrast to the global translation 
generator C[n], becomes second-class after introducing the new constraint B(z). 

With the Hamiltonian constraint, smeared out with a (lapse) function N, the constraint 
B has the Poisson bracket 



{B { z),H[N}) = N ( {K ' E "j E K » EV) B - 2vW) . 



(24) 



B will thus be preserved under the evolution generated by H, only if we add a new constraint 

A = 0. (25) 
The constraints A and B form a second-class conjugate pair 

{A(z),B(z')} = k5(z- z'). (26) 
The new constraint A must be preserved as well. It is diffeomorphism invariant in the 



full sense, and thus translation invariant, since 



{A(z),C[n]} = [n(z)A(z)]' 
is weakly equal to zero. The Poisson bracket of A with the Hamiltonian constraint is 



(27) 



{A,H[N]} =kN 

-kN 1 



OH 
~d£ 



dH 
~d£' 



dH 



'dH 


( m ) 




\d£" J 



kN" 



dH 
d£> r 



(28) 



The derivatives are 



dH 
~d£ 



H 
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+ - 



kVe 
i (e x e»)' 



(K X E X + K y E y )A 



1 / E x E y ' — E x 'E y 

E~ x E~y 



£ 



2 E x Ev 



£' - S" 



dH 



d£> k^E 



E x Ey 



and 



dH 



.£. 



dS" k^E 

With the new constraints A and B taken into account, we have the following weak equiva- 
lences 



dH _ E X E V " + E x "E y dH 
~d£ ~ ~KE x E~y ' dS 1 



0. 



Inserting this into ( 128|) gives 



iV 



(29) 



(30) 



So we see that, as in the case of B with diffeomorphisms, A is not invariant under the local 
action of H(z), so the full local Hamiltonian constraint becomes second-class, like the local 
diffeomorphism constraint C(z). If we choose a lapse function N linear in z and introduce 
the further constraint 

D := E x "E y + E x E y " (31) 

then the constraint A is preserved under evolution. 

For reasons that become more clear in later calculations, we make the more specialized 
choice iV = N(t). With the additional constraint d z N = on the Lagrange multiplier 
iV the system remains consistent, as can be checked using the constraint algebra. With 
this assumption the Hamiltonian constraint is reduced to a global condition H[N] = 0; 
the associated symmetry transformation is an evolution in a global time. This choice is in 
accordance with the form of the metric ([TJ, where the choice iV = 1 is even more special 



26] 



The constraint D is interesting. Imposing B = 0, we may express 



E x = Le~ p and E y = Le 13 



(32) 



and, after insertion into ( 131]) . the constraint equation D = becomes 2L(d].L+(d z f3)' 2 L) = 0; 
this is the spatial part of the Einstein equation ([2]). 
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Now D Poisson-commutes trivially with A and B. Its Poisson bracket with C[n], given 
n' = 0, is 

{D[f},C[n]} = -nD[f] (33) 

and vanishes weakly. However, the analysis is not complete since we have to be sure that 
D = is preserved under the Hamiltonian constraint. 

E. Consistency of the reduced system 

Taking into account the constraints A and B, and under the condition N = N(t), the 
Poisson bracket of D with H is, after integration by parts, 

{D[f),H[N}} « J dz f [N(K X E X " + K y E v " + K'^E* + K'^E y )] . (34) 

So far we have reduced a system on the six dimensional phase space with two first-class 
local constraints, H(z) and C(z), to one with five second-class local constraints H(z), C(z), 
A(z), B(z) and D(z) and two global evolution generators H[N(t)] and C[n(t)]. Numerically 
five constraints per space point would suffice to reduce six phase-space functions to one, 
corresponding to the free function j3 in ([2]), but consistency under time evolution requires 
more. 

Even with the assumption that N is independent of z, D does not weakly Poisson- 
commute with H[N]. The bracket is equivalent to 

- N J dz f(K"E x + KyE y + K X E X " + K y E y ") =: -N J dz f(z)J(z). (35) 

The new constraint J(z) can be expressed as a sum of similar terms including the derivative 
of D, 

J = E X [K X - E y ')" + E v (K y - E x ')" + E X "(K X - E yl ) + E y "(K y - E x ') + D' , (36) 
or, alternatively, 

J = E X [K X + E y 'f + E y (K y + E x ')" + E X "[K X + E yl ) + E y "[K y + E x ') - D' . (37) 
The Poisson bracket of J with H is weakly equal to the second derivative of D (using A and 
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B in the equivalence) plus additional terms 

{J[f],H[N]}n N J dzf^D" -2^ + ^ (E x 'E y/ — K x K y )' 



+2 



E x '\ 2 ( E y ' x 2 



E x J \Ey 



(38) 

(E x 'E yl - K x K y ) + A(E x "E y " - K' x K' y ) ' 



We clearly need to check the Poisson bracket of {J, H} with H. The constraints descended 
from J contain higher and higher derivatives so this leads to an infinite tower of constraints; 
the system in this form is inconsistent. On the other hand, we know from the metric 
formulation of (CQ) and (TSJ) that there is a consistent formulation for non-colliding waves with 
one configuration degree of freedom per point in light-cone coordinates. Obviously, there 
must be relations between the constraints to reduce the number of independent ones. 

This observation suggests an obvious solution to the apparent inconsistency. We can 
restrict the phase space variables at the kinematical level so that they only support left- or 

n 

right-moving waves [27j. The constraint J in the form ( |36l) or ( |37l) weakly vanishes when 
K x = ±E y> and K y = ±E X> . Then J and {J,H} are essentially D' and D" and so also 
weakly vanish. Hence we impose either the "right-moving" 

U x := K x - E y ' = and U y := K y — E x ' = (39) 

or the "left-moving" 

V x := K x + E y/ = and V y := K y + E x ' = (40) 

as primary constraints. As shown in Appendix A these relations, together with the equation 
of motion, provide a consistent solution to the Einstein equations in terms of the triad 
and canonical momenta. The relations also anticipate that the canonical momentum of the 
metric variable E p is equal to ± its spatial derivative for pp-waves. 

In the following we work with the right-moving constraints U p . The Poisson brackets of 
U x and U y are 

{U x (z),U y (z')} = {U y (z), U x (z')} = 2k5'(z - z') = 2k^- 5{z - z') 

oz (41) 

= -{U y (z'),U x (z)}. 

Note the antisymmetry in z and z' in spite of the symmetry under the exchange of U x and 

Uy\ 
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These right moving constraints have non-vanishing Poisson brackets with B 

{U x (z), B(z')} = KE y (z)5(z - z'), and {U y (z), B{z')} = kE x (z)5(z - z'). (42) 

Introducing the multipliers u p (not to be confused with the light-cone coordinate u in the 
introduction) for the constraints U p and h for D we have the Poisson bracket 

{U x [u x ],D[h]} = K J dz{[u x (z)h(z)\" + u x (z)"h(z)}Ey(z) } (43) 

and a similar relation for {U y [u y ], D[h]}. The Poisson brackets {U p , H[N(t)}} and 
{Up, C[n(t)}} vanish weakly. Thus, the constraints U p are compatible with time evolution 
and their introduction solves the problem of the infinitely many constraints, thus making 
time evolution consistent. On the other hand, this introduction increases the number of 
second-class constraints to seven, which is definitely too many. What remains to solve is 
this apparent overconstraining of the system. 

Physically the reason for the constraints U or V lies in the fact that the full Hamiltonian 
constraint of plane gravitational waves applies to modes going both into the positive and the 
negative z direction and their mutual interaction. A superposition of left- and right- moving 
waves would introduce complicated interactions and spoil the simple form of the metric. 
In section IIIIDI we will see that under the conditions U p = or V p = the Hamiltonian 
constraint generates simple plane wave propagation. 



III. DIRAC BRACKETS 



In this section we construct the Dirac brackets of the local second-class constraints step 
by step (see below), according to algebraic relationships. The algebraic Poisson bracket 
structure associates the second-class constraints into two "pairs", (A,B) and (U x ,U y ), and 
three single constraints D, H, and C. In addition A and B do not contain derivatives 
and so are actually associated to each point z separately. In the course of the analysis the 
constraints C and H turn out to be dependent, more precisely, equivalent to D, so that the 
set of independent constraints reduces to the convenient number of five. An odd number 
of second-class constraints (per space point z) may appear incompatible with the standard 



construction of Dirac brackets 



211 ]. but not all of them are related exactly to one point, 



some of them contain derivatives. 
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For a mechanical system with second-class constraints Cj, i = 1,2, ...,2n the Dirac 
bracket of two phase space functions F and G is defined as 

{F, G} D = {F, G} - {F, Ci}M£{C k , G} (44) 

in terms of Poisson brackets. The matrix M^} is the inverse of the matrix = {Cj, Ck] 
of the Poisson brackets among the constraints. After the Dirac brackets are constructed, 
the constraints can be imposed strongly. This reduces the system to its physical degrees of 
freedom. 

A helpful fact about Dirac brackets is that they can be constructed recursively, i. e. the 
construction of equation fl44P can be carried out for any subset of second-class constraints, 



provided the matrix of their Poisson brackets is invertible [2l|]. After imposing these con- 
straints strongly the procedure can be repeated with the preliminary Dirac brackets replacing 
the Poisson brackets in equation ()44l) . This possibility greatly facilitates the work with our 
constraints. In field theory, of course, the matrix multiplication in (jHJ) implies integration. 



A. Dirac brackets, version Dl 

Beginning with the pair (.4, B) we have the Poisson brackets (I2b'p and 

Mr k \ z ,z') = -5(z-z')\ 01 I. (45) 

The ensuing Dirac brackets, version Dl, are explicitly 

{F(z),G(z')} m = {F(z),G(z')}-^ J dz"{F(z),A(z")}{B(z"),G(z>)} 
+i J dz"{F(z),B(z")}{A(z"),G(z% 



(46) 



Due to the appearance of A in both the integrals on the right hand side, neither the Dirac 
brackets of the variables E p , K p , nor those of the remaining constraints differ from the 
corresponding Poisson brackets. We can simply impose A and B strongly. When this is 
done, U x , U y , and D are untouched, whereas C and H are simplified considerably: The 
diffeomorphism constraint drops its term — -£'A and becomes 

C = - (K' X E X + K' y E y ) , (47) 
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whereas the Hamiltonian constraint boils down to 

H = —— [K x K y + E x 'E y ' - (E x E y )")} (48) 

with the second-derivative term not contributing to integrals with a ^-independent test 
function. Without this term the last expression for H is similar to the Hamiltonian of 
two free Klein-Gordon fields, the non-linearity of GR is now hidden in D, which is not 
conserved under the evolution generated by the Hamiltonian constraint. In the simplest 
case of constant lapse and shift, e. g. N = n = 1, D commutes weakly with C, 

{D,C[1]} = E x E yl " + E x 'E y " + E x "E y ' + E x/ "E y = D', (49) 

but not with H, 

{D, H[l}} = E x Kl + E X "K X + E y Ky + E y "K y , (50) 
as long as we do not introduce the constraints U p . 



B. Dirac brackets, version D2 

The next pair of second-class constraints, (U x , U y ), has the mutual Poisson brackets (|4ip . 
To construct the inverse of the matrix of these brackets we need the inverse of the derivative 
of a (5-function, denoted by S^ 1 ^, which satisfies the relation 

J dz"5'(z - z")5^ l \z" - z') = 8{z - z'). (51) 

Obviously 5^~ x \z — z') is a step function plus an additive constant that is adjusted by 
demanding antisymmetry 28] 



<J (_1) (z -z') = - sign(^ - z') . (52) 



We construct the matrix N ik that plays an analogous role as M ik in 



N7 k \ z ,z') = ±sign(z-z') ( ° * | . (53) 
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With this matrix the next version of Dirac brackets 
{F, G} D2 = {F, G} - {F, U x }N; k l {U y , G}, becomes 

{F(z),G(z')} D2 = {F(z),G(z')} 

~ [ + dz"dz"'{F(z), U x {z")} sign(z" - z"'){U y (z"'), G(z')} (54) 

~ [ + dz"dz"'{F(z), U y {z")} sxgn{z" - z"'){U x {z>"), G(z% 

(The Dl brackets are the same as the Poisson brackets so the label is omitted.) In particular, 
the Dirac brackets of the remaining fundamental variables are the following 

{K x (z),K x (z')} D2 = {K y (z),K y (z')} D2 = 0, 

{K x (z),K y (z')} m = ~5'(z-z'), 

{K x (z), E*(z')} m = {K v (z), Ey(z')} m = %S(z- z') 

{K x (z),Ey(z')} m = {K y (z),E*(z')} m = 0, 

{E%z),E*(z')} B2 = {Ey(z),Ey(z')} m = 0, 

{E*(z),Ey(z')} m = ^sign(z - z'). 

The bracket relations between E x and E y may look awkward due to non-locality. This is 
explained by the form of the constraints U x and U y . Integrating them yields the E's in form 
of an integral over K. The expression 

1 



(55) 



E y (z) 
v 1 2 



r K x { z ')dz' - r K^dz' = i r ^ z - ^ 

*J Z 1/ Z J */ z 



and its counterpart E x from U y = make the non-locality of their Dirac brackets plausible. 
These brackets are of the same form as those of the self-dual fields of j^J, see section V. 

To impose the U constraints strongly, we can express the i^'s in terms of the _E"s or vice 
versa, or choose one of the canonical pairs (K x , E x ) and (K y , E y ) as fundamental variables. 
To preserve the canonical structure, the latter choices would seem to be preferable, but in 
different calculations different choices may be suitable. 

After the C/'s are imposed strongly, C(z) and H(z) become equivalent to -D(z), explicitly 

C = - D + - (U' X E X + U' y E y ) (56) 
16 



and 

H=-D-- (U x U y + U X E X ' + U y E»'), (57) 

K K 

this means that finally the number of independent local constraints is reduced to five and 
that D(z) now implies also the global constraints H[N] and C[n]. So the constraints U p 
themselves solve the problem of overconstraining that arose after their introduction. Further, 
the fact that U p (V p ) lead to H[l] = ±-D[l] confirms that U/V single out left/right moving 
wave modes. The integrated Hamiltonian constraint with N = 1 becomes 

H[l] = ~J dzE x \z)E y \z) = ~\J dzK x (z)K y (z). (58) 

Finally D commutes with the total Hamiltonian, which is now (for N = 1 and n = 0, 
according to the assumption in (CQ)) just H[l], 

{D(z),H[l]} D2 = D'(z)^0. (59) 
C. The final Dirac brackets 

At this point we have one phase space degree of freedom, represented equivalently by one 
of the above-mentioned pairs of variables, one local constraint D(z) per point z and one 
global one, H[l], which is at the same time the generator of time evolution. The constraints 
D(z) are second-class and their Dirac brackets, version D2, are rather complicated. 

{D(z), D(z')} D2 = k [f(z, z')S"'(z - z') + g(z, z')5\z - z') + h(z, z'^iz - z')] , (60) 

where 5^ l \z — z') was introduced in (|52p and 

f(z,z>) = l - [E%z)Ey(z') + E*{z')Ey{z)\ , 

g{z,z') = \[E*{z)Ey"{z>)+E*"{z)Ey{z') 

2 (61) 

+E x {z')E y "{z) + E x "{z')E y {z)\ , 

h(z,z') = X - [E x,, (z)E y "(z') + E x "(z')E y "(z)] . 

Let's denote by A(z,z') the inverse of {D(z), D(z')} D2 , needed in the construction of the 
final Dirac brackets 

{F(z), G(z')} D = {F(z), G(z')} m - (62) 
(k)- 1 J dz" dz"'{F(z), D(z")} m A(z", z'") {D(z"% G(z')} m . 
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It is an antisymmetric function satisfying 

dz"{D{z),D{z")} D2 A{z",z') =5{z-z'). (63) 

We do not have the full solution to this equation. In Appendix [B] we calculate an approxima- 
tion, demonstrating some qualitative features of the canonical structure rather than giving 
the exact Dirac brackets. In the following A is understood as this approximation and { , }d 
as a representative part of the full Dirac bracket, constructed with A. 

As already mentioned, when we apply five local constraints strongly, there remains one 
free variable. If we choose E x , our fundamental Dirac brackets are those of E x at different 
points, constructed according to (I62I) . For this purpose we need the bracket 

{E*(z), D(z")} m = E*(z") {E*(z), V')>D2 + E x " (z"){E x (z), ^(/)} D2 

(64) 



I E x (z") S'(z - z") + I E x "{z") sign(z - z")- 



{D(z"'), E x (z')}d2 is calculated analogously. For our approximation of {E x (z) , E x (z')}b we 
take only the 5' parts of these brackets. In the following A, calculated in fIBlip . is more 
conveniently expressed in terms of anti-derivatives of the 5 functions, 

8(~ 3 \z - z') = \\z - z'\(z - z'), 6^(z - z ') = ±\z- z'\(z - z'f. 



Putting these ingredients together, we have 



{E x (z), E x {z')} D « - / dz" dz" 1 E x {z") 5'{z - z") 



5(-z)(z"-z"') 
A(^) 



ii \ -,1" \ \ 2 



3(A(^)) 



8^\z" -z'") 



E x (z w )S'(z'" -z'). 



4 

After integrating the 5' functions by parts we expand E x (z) and its first derivative around 
z = (z + z')/2 



E x (z) « E x (z) + E x '(z 



z — z 



+ ... 



and E x (z'), and make use of 



S^\z - z') ■ (z - z') 2 = 5^ 2 \z - z') ■ (z - z') = 2 5^\z - z'). 

With all the variables evaluated at z (so that A corresponds to A in Appendix B) and 
inserting finally A = L 2 and E x = Le~P, we find 



{E x (z), E x (z')} D « - e~ 2p sign(^ - z' 



l-[2^') 2 + ^")(z-^ 2 



(65) 
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In the flat space-time regions z\ < z < Z- and z + < z < z^ of our coordinate domain, where 
(3 = (and L" = 0) the field E x satisfies bracket relations analogous to the commutation 



relations of self-dual Klein-Gordon fields, considered in [22J, which are constructed by re- 
striction to waves going into one direction. The correction in the brackets for E x is expressed 
purely in terms of the wave factor /3. Although the approximation is rather qualitative, it is 
quite instructive for some insight into the influence of the gravitational Hamiltonian in the 



canonical structure of the self-dual fields 



22 1 



For (3 the above expression can be interpreted as a low-order approximation of 
a gravitational correction. Were this quantized, this would appear as a variable Planck 



constant, as suggested by Hossenfelder [23], or a variable gravitational constant. Other 
corrections, however, do not fit into this scheme, they give rise to qualitatively different 
terms from f )65|) . 

D. Time evolution 

The time evolution of a phase space function F is generated by its Dirac bracket with the 
total Hamiltonian. As already stated in the preceding section, by virtue of the U's H[N] 
becomes equivalent to C[N] and so the total Hamiltonian becomes a generator of a rigid 
translation. (Had we chosen the V constraints, H[N] would be equivalent to —C[N].) This 
equivalence allows to introduce C[N] as a true total Hamiltonian, when we choose a lapse 
function. The most convenient choice N = 1 means a constant unit of time. 

The Hamiltonian constraint H tot = H[l] = C[l] being first-class at every stage, its Dirac 
brackets with any phase space function are equal to the corresponding Poisson brackets, 
which are equal to the ^-derivative, according to the nature of C[l] as translation generator. 
Hence, 

F(x) = {F(z), H[1]} D = {F(z), C[l]} = F\z). (66) 
Equivalence of H[N] with ±C[iV] simply means that time evolution is a rigid space trans- 



22] 



lation to the left or to the right, the same relation that characterizes self-dual fields 

This completes the Einstein equation by making all the variables depend on t — z (or 
t + z, alternatively). So we have recovered the classical equation of motion fl2]) in a much 
reduced phase space. One can describe the system with a single function, e. g. E x (z), on 
(zi,Z2) subject to the second-class constraint D. 
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IV. THE IMMIRZI PARAMETER AND THE POLARIZATION ANGLE 



In this section we return to the orthogonality of the connection components ( Tl8|) to show 
that this is satisfied by the reduced model. After strong imposition of all the constraints 
the two-vectors A x = (A^ ,A X 2 ) and A y = (A y 1 ,A y 2 ) are orthogonal and there arises a 
simple relation between the angle a between E x = (Ef, E%) and A x = (A x , A 2 ) and the 
corresponding angle a. 



The variables in polar coordinates from 16j, with (3 = (3 = 0, corresponding to the gauge 
£ = i] = are 

Ef = E x , El = 0, E\ = 0, E\ = EV; 

(67) 

A\ = A x cos a, A\ = A x sin a, A y = —A y sin a, A y = A y cos a. 

From elementary calculations of the connection components T in terms of the E's we find 
in the gauge j5 — /3 — 

r^ = r 2 = o, rl = -E»' t t\ = -e*\ (68) 

so that the diagonal components A l x and A y contain only extrinsic curvature, 

A* = A x cos a — jKl, A 2 y = A y cosa = ^K*. (69) 

On the other hand, from the Gaufi and the polarization constraint {£, H} = we obtain 

(H, (A.14)) 

A\ = Aj-sina = T x , A y = —Ay sin a = — T y , (70) 

thus the off-diagonal components are purely composed from T's. Now the vectors A have 
acquired the form 

A x = (jK x , -&'), A y = (-E x ',-fK y ), (71) 

where we have written K x = K x and K y = K y . For the absolute squares of these vectors 
we get 

(A x ) 2 = ^ 2 (K x ) 2 + (Ey') 2 and (A y ) 2 = ^(K,) 2 + (E x ') 2 . (72) 
With the constraints U p (or V p ) this becomes 



A x = K x y/T+i* and A y = K y ^/l + 1 2 . (73) 
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Inserting into ([67)1 and comparing with ( |69|) . 

A x x = K x ^l + j 2 cos a = jK x , A 2 y = K y y/l + i 1 cosa = jK y , 

leads to a relation between the angles a and at and the Immirzi parameter: 

a = a = arccot7, (74) 

so in the end, after all gauge fixing, A x and A y are orthogonal and orthogonality is compatible 
with the Hamiltonian constraint, when the latter reduces to a translation generator. 

V. SUMMARY AND CONCLUSION 

The principal aim of our considerations is the loop quantization of polarized gravitational 
plane waves to see whether dispersion relations would be modified, or if there are other 
effects from the granularity of the kinematic states such as a variable speed of gravitation 
or a variable Planck constant. One way to handle a quantum theory of plane waves is to 
quantize a more general system, such as a model with plane symmetry, analogous to the 
Gowdy model exploited here, and then to distinguish the subspace of left- or right-going 
modes of the full Hilbert space. But, as the formulation of basic operators in |l7J shows, 
this turns out to be quite complicated. 

In the present work we reduced the formalism of plane waves to the physical degree of 
freedom at the classical level. To derive a classical description of plane waves suitable for 
loop quantization we started with the assumption of homogeneity of the spatial geometry 
in the transversal directions and a coordinate system extending in both directions beyond 
a gravitational pulse, so that the latter one is embedded between two slabs of flat space 
in these coordinates. The finite range of this coordinate system rendered the integrations 
over the remaining spatial coordinate finite. This setting is analogous to the Gowdy model 
described in |16[ and we used the formalism developed therein. The reduction was completed 
with a Dirac constraint analysis. 

This was done as follows. The description starts with a symmetry reduced model with 
three configuration space degrees of freedom and the standard diffeomorphism and Hamil- 
tonian constraints of GR, giving one field degree of freedom per spatial point. The first step 
of pp-wave reduction was carried out in section III Dt We restricted the metric in the triad 
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variables to the simple diagonal form of (JTJ) by introducing the constraint B. Preservation 
of this constraint led us to introduce the secondary constraints ^4 = and D = 0, the latter 
one b e lng the (spatia! p roJ ect 10 n „ t ) the classica! E inStei n e q na tl on Q. Knewin g that the 
symmetry reduced system cannot describe colliding waves and requiring consistency in the 
sense that further secondary constraints vanish, we introduced the right- (or left-) moving 
constraints U p (or V p ) in section III El After restricting the fields to left- or right-moving 
modes by imposing constraints we find a system analogous to the self-dual fields - also 
scalar fields propagating in one direction - described by Floreanini and Jackiw 221 ]. 

We impose these constraints in a special coordinate system, fixed by the constraint B, 
which breaks diffeomorphism invariance, but makes the form of the constraints U p simple. 
Alternatively, the constraints distinguishing left/right going modes can be formulated in a 



25]. 



fully diffeomorphism invariant manner. Work is underway on this approach 

In the second step of the reduction we used the "pair-wise" structure of the Poisson 
algebra of the constraints to recursively construct Dirac brackets. In section HTT1 using A and 



B and then U p we constructed the first two versions of Dirac brackets. After imposing the 
constraints U p , the constraint D(z) commutes with the Hamiltonian, but becomes second- 
class by itself. At this point the evolution has become trivial; all the complications are now 
in the D2-bracket relation of D(z) and D(z'). To complete the canonical treatment of the 
problem, we constructed the final Dirac brackets. In the end the D(z) are the remaining 
second-class constraints, leading to one variable (we have chosen E x ) with a very non-trivial 
Dirac bracket {E x (z), E x (z')}d, containing the step function sign(,z — z'), multiplied by an 
analytic function in z — z' . An approximation is given in equation (1651) . It is clear from 
this bracket that the canonical structure of the reduced system is obscure when using this 
variable. 

This makes a full quantization, i. e. a formulation in terms of operators and a Hilbert 
space, elusive - we do not even know the exact closed form of the Dirac brackets. The 
lowest-order terms, however, turn out to be the analog of the commutation relation of a 
linear self-dual field 22] plus gravitational corrections. 

The result of our preliminary classical considerations gives no suggestion of dispersion in 
these waves, which would provide an indication of an energy- dependent speed of gravitation. 
The reduction to left or right moving waves leads automatically to the equivalence of the 
Hamiltonian to the generator of spatial translations. This equivalence was not assumed from 
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the beginning, as in other approaches, for example in light-like coordinates, but it appears 
as a result of the disentangling of otherwise colliding modes by the constraints U p . In this 
way this complete reduction with the aid of Dirac brackets differs from other canonical 
approaches such as [24]. In a quantum theory, based on our classical analysis, where these 
generators are promoted to operators, an analogous result can be expected. The final Dirac 
bracket does hint at a modification of the quantum relations. 

The non-local bracket of E x (z) fl65l) suggests a modification of the Planck constant (or the 
gravitational constant) in the first approximation, rather than a variable speed of light. In 
the framework of our approach (starting from unmodified GR) we expect that the space-time 
texture arising in a quantum theory of gravity would influence the fundamental structure of 
quantum theory, mainly the commutators and the uncertainty relations derived from them. 

Finally, a remark on the approach to the reduction: After specializing to a z-independent 
shift vector and a lapse function of the same type in order to conserve B and the ensuing 
constraints we could have abandoned the local constraints C(z) and H(z) and kept only the 
according global ones C[l] and H[l]. This would mean to start with a theory different from 
GR, a theory without full diffeomorphism and time-reparametrization invariance. Neverthe- 
less, we would have arrived at the same results, because on the constraint surface determined 
by A, B, U x and U y the local constraints H(z) and C{z) are equivalent to D(z). For this 
reason we did not have to make explicit use of the local diffeomorphism and Hamiltonian 
constraint in our work. 

Appendix A: The Ricci tensor for the metric diag (—1, (E y ) 2 , (E x ) 2 , l) 

To introduce canonical variables, we calculate from the Levi-Civita connection the ex- 
trinsic curvature components 



as canonical conjugate variables to the metric variables E x and E v . In terms of these 
canonical variables we get the five non- vanishing independent components of the Ricci tensor: 



k x = Ey, 



K, 



y 



E x . 
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Ru = E y (K x - E y ") + —(K x K y - E x 'E y '), 

R 22 = E x {K y - E x ") + ^(K x K y - E x 'E yl ). 

The vacuum Einstein equations i?03 — and -R33 = are the constraints C and D, divided 
by —E x E y , the remaining ones contain time evolution. Here C and D appear as primary 
second-class constraints, unless we smear out C with a ^-independent function. As in the 
main text, all the constraints are consistent when U p = or V p = and d/dt = ±d/dz. For 
a quantization, also this approach requires Dirac brackets. 



Appendix B: The inverse of {D(z), D(z')}d2 



In the defining equation (|63|) of A(z, z'), with fl64|) inserted, the derivative of the 5 function 
is shifted to its second argument: 8'{z — z") = —d/dz"5(z — z"), and analogously the third 
derivative. Then, integrating by parts, we get 



dz" {£^3 [/(*, *')\ + ^ \9(z, z")A(z", z')\ \ 5(z - z") 

+ f dz" h(z, z") 5 { - 1 \z - z") A(z", z') = 5(z - z'). 



(Bl) 



Consider the third derivative of the first square bracket, 



d 3 \ ( d 2 \ d 



d \ d 2 d 3 



dz" 3 J \dz" 2 J dz" \dz" J Jdz" 2 dz" 3 

This is multiplied by 5(z — z") in (IBip . so we need / and its derivatives at z" = z. 



(B2) 



f(z,z) = E x (z)E y (z), 



df 



d 2 f 



dz 



1/2 



dz" 
d 3 f 



(E x (z)E y (z)Y. 



dz 



1/3 



= - (3D'(z) - (E x E y )"') . 

-z 

= ^D'(z), so the second square bracket 
does not contribute anything, when D w 0. Denoting E x (z)E y (z) = L 2 (z) by A(z), we may 
write (IBlj) in the form 



Similarly g(z,z) = D(z) and d/dz" g(z, z") 



d 3 3d 2 1 

A(z) — A(z, z') + - A'(z)—A(z, z') - - A'"(z)A(z, z') 

+ / dz" h(z, z'^-^iz - z") A(z", z') = 5(z - z'). 



(B3) 
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In this equation A can be considered as a Green function of an integro-differential oper- 
ator. Because of the required antisymmetry of the Dirac bracket we are looking for an 
antisymmetric Green function. 

The dominant coefficient function is A, the square of the background factor of the grav- 
itational wave, its derivatives are smaller. To find an approximative part of the solution of 
this equation we first look for a solution of the differential part and leave the integral part 
for later iterative corrections. The leading term contains a third derivative, therefore the 
leading term in A is expected to contain the function 



the third derivative of which is the delta function on the right hand side. To find at least 
an approximation for A, we make an ansatz in the form of a product of this function by 
a symmetric function of z and z' . Further we assume this function to be analytic in some 
neighborhood of z — z' . Written in terms of z + z' and z — z', this function has only even 
powers in z — z', so that 




(B4) 




(B5) 



with the coefficients 




For the derivatives with 



respect to z we find 




A, zz (2, z') 



sign(^ — z') 



4 



2a 2 + 2a' 2 (z - z') + 




+ 12a 4 - z'f + ... 



A,zzz(z, z') = a 2 5(z - z') + 



sign(z — z') 



4 





Inserting this into the first line of ( IB 3 1) gives 



Aa 2 5(z - z') + -(Aa 2 + A'a 2 ) sign(z - z') 




(B6) 



+ (^Aa' 2 ''+J-AV 2 ' 




' 2 ' - 1 A'"a 2 + ^ Aa' 4 + ^ A'a 4 ] \z - z'\(z - z') + . . . 
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Note that A and its derivatives are functions of z, whereas the a» are functions of z and their 
derivatives refer to this argument. 

To express everything in terms of z and z — z', we expand A(z) = K(^f- + ^^-) around 

z: 

A(z) = A + ~ K(z - z') + ~ K(z - z'f + ... (B7) 



where A = A(z). Inserting into (1B6|) and rearranging terms gives 



3 

A a 2 5(z - z') + -(A a 2 + Ao<2 2 )sign(,2 - z) + 



/3 9 3 \ 

( g A a 2 + - A(,a 2 + - , Ag'a 2 + 6A a 4 j \z - z'\ + 

l - A a 2 " + ^ A' a 2 ' + ^ A' 'a 2 + A' "a 2 + ^ A a 4 + y A[,a 4 ^ x 
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\z-z'\(z-z') + ... 

The coefficient of the 5 function on the right hand side of ( 1B3I) is one, therefore Aoa 2 = 
1. Now all primes denote derivatives with respect to z, so we can split off the vanishing 
derivatives of A a 2 and have 





/I 




- z 



m 

/3 3 9 15 \ 

( g A 0«2 + 1 + 2 A ° fl 4 + y A a 4 ) I* " A{* -/)+••• 

Inserting now a 2 = Aq 1 and setting the coefficient of |z — z'\ equal to zero we find 

a ~ 1 (A ' o)2 fB9^ 

° 4 - i6^r (B9) 

or, in terms of L, 

- = lk' °* = i W (B10) 

The coefficient of |z — z'\ {z — z') in the last term of ( 1B8I) cancels, so our ansatz ( 1B5I) leads to 
an antisymmetric approximation of the inversion of {D(z), D(z')}^ 2 in some neighborhood 

of z = z' 



A(z,z') « 



\ + l -^±(z-z<y 



(Bll) 



4L 2 (z) 

This is an approximation to the differential part of the integro-differential equation ( 1B3I ). 
iterating the integral part would not give contributions of this order, but only higher anti- 
derivatives of 5. 
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On the other hand, had we begun with the integral part, we would have obtained 8' as 
first approximation, because the integral contains the first ant i- derivative. Beginning with 
this, we would get higher and higher derivatives of the 5 function in the sequel. This does 
not show up in our A and so the ensuing expression for the Dirac brackets reveals only part, 
although an important part, of the consequences of imposing D strongly. 
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The choice N = N(t) is also essential for calculations of the algebra of constraints including 
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equation (j3l|) and {U p , H[N]} « 0. The constraints C/ p are defined in equation (f39j) . 

There is precedent in the literature for this. Scalar fields in 1+1 dimensions with the same 

property are called "self-dual fields" and are quantized, for example, by Floreanini and Jackiw 



221 ]. The restriction of gravitational field variables to either left or right moving waves will 
lead to analogous bracket relations as those for the self-dual fields with corrections in the part 
of space-time where the wave factor /3 is non-constant. 
Up to boundary terms that do not apply when z is in the interval (zi, z%). 
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